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A SATELLITE ALTIMETER BIAS RECOVERY SIMULATION 


ABSTRACT 

Simulations are made for bias recovery of a satellite- 
borne altimeter by means of data from various combina- 
tions of ground-based range and angle measuring systems 
to determine the accuracy with which the altimeter bias 
can be recovered under these restricted conditions. Among 
the parameters considered are the number of stations, 
station geometry, elevation angle of the pass, and the re- 
lative amount of range or range with angle data used. An 
assessment of the probability of obtaining adequate laser 
and camera data due to the effect of weather and day or 
night conditions is included. 

The results indicate that, for the systems considered 
in these simulations, the two station network, with collocated 
lasers and cameras, is best with respect to accuracy, 
operational data collection probability, and minimum cost. 
Such a network is therefore recommended for the satellite 
altimeter calibration (GEOS-C). 
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A SATELLITE ALTIMETER BIAS RECOVERY SIMULATION 


I. PURPOSE 

This altimeter simulation study was initiated to determine how accurately 
the altitude of a satellite could be determined over short arcs with data from a 
satellite -borne altimeter and from various combinations of ground based range 
and angle tracking systems. The tracking stations were chosen to take advan- 
tage of existing sites or sites that could be readily equipped. The us. of ranging 
data was assumed at the outset to be essential. Thus, in this simulation, there 
are never more angle trackers than range trackers. Also the angle trackers, 
when used , are always collocated at one or more of the ranging sites , as would 
most likely be the case if the technique is implemented. 

The range data a priori bias and rms noise values are chosen to conserva- 
tively represent values which were obtained from several lasers and C-band 
radars during the GEOS-2 tracking system intercomparisons (references 1,2, 

3, 4, 5, 6). The angle data a priori bias and rms noise values represent, some- 
what optimistically, values obtained from several cameras intercompared on 
GEOS-1 (reference 7). 

When a Goddard laser was collocated and intercompared with two C-band 
radars at Wallops Island, (references 1, 2, 3) the results indicated that those 
C-band radars, when certain calibration and data preprocessing errors were 
removed, were more or less comparable to that laser in accuracy. Over the 
last year, however, the laser data rms noise and accuracy has improved by a 
factor of 2 or 3, and over the next few years is expected to further improve. The 
C-band radar does have the operational advantage of being able to observe through 
the clouds. However, in the laser and C-band intercomparisons, the occasional 
large discrepancies which arose between the systems were always accounted for 
by some problem with the C-band and not the laser. 

From this experience, it seems advisable for an altimeter calibration op- 
eration to utilize both the C-band and the laser, the C-band for its operational 
advantages and the laser to verify the C-band results. For convenience, in this 
report, the range and angle systems are referred to as lasers and cameras. 

This simulation study attempts to answer the following questions: 

1. How is the accuracy of the altimeter bias recovery affected by number 
of lasers and cameras tracking ? 

2. Is it necessary to have camera data or will the less accurate laser angles 
provide the directional information ? 
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3. How is the altimeter bias recover}’ affected by the elevation angle of 
the satellite pass ? 

4. How sensitive is the altimeter bias to location of the laser stations? 

It is recognized that the answers derived are highly dependent on some of 
the a priori assumptions, such as the assumed 2 meter ranging accuracy (bias) 
for the laser. However, the a priori values chosen were thought to be fairly 
realistic, based on the tests described in the references quoted on page 1. 

More recent developments indicate that a laser range accuracy of 0. 5 meters 
might now be more realistic. 

This study is not intended to provide comprehensive answers to the above 
questions , but it is hoped the results will be useful in helping to establish a plan 
for evaluating a satellite -borne altimeter having an overall systems accuracy 
requirement of 5 meters or better, (reference 8). 

Part of the 5 meter overall systems error budget is used up in external er- 
rors such as a refraction error, an error due to sea state, and an error due to 
non-nadir reflection from the sea surface. Consideration of these errors means 
the calibration system must determine reference heights of the spacecraft alti- 
meter above mean sea level (MSL) to better than 5 meters. For the assumptions 
made in Reference 9, the required accuracy is 4. 1 meters. 


II. APPROACH 

The study is divided into four cases to answer the four specific questions. 

In general , the inputs to the computer program are kept either as unconstrained 
or as realistic as possible. These inputs are such quantities as the orbital ele- 
ments, station locations, observation uncertainties , observation sampling rate, 
and error model uncertainties. In all cases, the same orbital elements and ob- 
servation uncertainties are used. 

The four cases are as follows: 

Case 1 : 

The number of laser stations is varied from a minimum of one to a maxi- 
mum of four. For each number of laser sites, the number of cameras is 
varied from zero to the number of lasers. The number of cameras is al- 
ways less than or equal to the number of lasers. The situations are sum- 
marized in the following matrix: 
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CAMERAS: 


TRD 

CUR CUR 





GTK 

GTK 

GTK 



ANT 

ANT 

ANT 

ANT 


0 

1 

2 

3 

4 

LASERS: ANT 1 

(X) 

(X) 




GTK, ANT 2 

X 

(X) 

X 



CUR, GTK, ANT 3 

(X) 

X 

X 

X 


TRD, CUR, GTK, ANT 4 

X 

X 

X 

X 

(X) 


For each situation , recovery is made for the orbital elements ; a constrained 
range bias for each laser, the constrained survey for each station, and the 
altimeter bias. 

Case 2 : 

The number of laser stations is varied from one to four. Tracking is sim- 
ulated with laser azimuth and elevation angles added to the laser ranges in 
order to determine what advantage the addition of the laser angle data pro- 
vides over the range data alone. 

Case 3 : 

a. Tracking is simulated for one laser and one camera (1L1C) collocated at 
the Antigua station. For this station the satellite pass rises to a maximum 
elevation angle of 87° . Above 20° , the pass is broken into 14 one-minute 
mini-arc spans of data. Each mini-arc is considered as a separate pass 
containing both laser and camera data , and recovery is made for the orbital 
elements, a constrained laser range bias, the constrained station survey, 
and the altimeter bias. This case, in effect, simulates a geometric recov- 
ery for the altimeter bias for various elevation angles, since the mini -arcs 
contribute very little dynamic constraint to the solution. 

b. An altimeter bias is recovered for most of the 14 mini -arcs as in Case 
2a. However, the orbital elements, laser range biases, and station surveys 
are recovered for the total arc. This case shows how the altimeter bias re- 
covery depends on elevation angle when only short spans of altimeter data, 
but longer spans of comparison data are available. The station geometry is 
varied to include the minimum 1LOC and the maximum 4L4C situations, the 
1L1C situation as in Case 2a, and the intermediate 2L1C and 3LOC cases for 
comparison as indicated by the () symbols in the matrix given in Case 1. 
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c. The error in altitude for each of the mini -arc epoch points used in Case 
3a is calculated by means of a geometric error propagation using range and 
angle data from ANT. This error propagation assumes a 2m range meas- 
urement uncertainty, a variety of angle measurement uncertainties, and the 
same survey uncertainties as in Case 3a. 

Case 4 : 

a. A 3-laser, 0-camera (3LOC) situation fora widely spaced laser network 
having laser stations at Antigua, Key West, and Panama, as was used in Ref- 
erence 10, is compared with the3LOC situation in Casel having more closely 
spaced lasers at Antigua, Grand Turk, and Curacao. 

b. A 2-laser, 0-camera (2LOC) situation with lasers at Antigua and Trinidad 
is explored using more optimistic horizontal survey and laser range noise 
and bias constraints and a higher laser data rate. This is the only case 
where these survey and laser data constraint parameters are changed. 


III. CONSTRAINTS 

Realistic station locations and orbital elements are used and the a priori 
estimates of the uncertainties in the orbit, the survey, the measurements, and 
the system error model terms, are generally either realistic or conservative. 
The same values are used for all cases except case 4b. The simulated data are 
generated at realistic intervals for both the lasers and the cameras. The par- 
ticulars are as follows: 


1. STATIONS 


The stations are chosen in the Caribbean Sea area where tracking sites al 
ready exist or where they could be easily installed. 



Latitude North 

Longitude East 

Height (MSL) 
(Meters) 

Antigua 

17° 08» 38. 5” 

298° 12* 25.7" 

58 

Grand Turk 

21° 27 f 46. 4" 

288° 52’ 4.4" 

28 

Trinidad 

10° 44 f 35. 5" 

298° 23 T 28.9" 

792 

Curacao 

12° 05’ 26.2” 

291° 07* 43.9" 

5 
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The following two stations are also used in Case 4a: 


Key West 

24° 0 f 

0” 

278° 

0* 

0" 

0 

Panama 

CD 

o 

© 

0*' 

281° 

O’ 

0" 

0 


The stations are ordered for the several cases as follows: 

1 Station: Antigua (ANT) 

2 Stations: Antigua, Grand Turk (GTK) 

3 Stations: Antigua, Grand Turk, Curacao (CUR) 

4 Stations: Antigua, Grand Turk, Curacao, Trinidad (TRD) 


In the instances where there are both cameras and lasers at the same *lte, 
they are assumed collocated at the same position. 


2. ORBIT 


The orbit assumed is representative of the type of orbit originally consid- 
ered for the GEOS-C satellite. However, any trajectory providing a similar 
tracking geometry will produce similar results. 

Osculating elements: 


Semi -major axis, a (meters) = 

Eccentricity, e 

Orbital plane inclination i, (deg) = 
Argument of perigee, cj (deg) = 

R.A. of ascending node, £2 , (deg) = 
Mean anomaly, M, (deg) = 


7817906.7 
0.025 
19.996 
0.251 
-0. 033 
53. 844 


The orbital elements and their a priori constraints in earth-centered, earth- 
fixed Cartesian coordinates are: 


EPOCH 

E(m) 


Elements Constraints 

= 35460 seconds after 0 h GMT 

= 3523131.728 1000 meters 
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Elements 


Constraints 


= -6489174.888 


F(m) 
G(m) 
E(m/sec) 

F(m/sec) 

• 

G(m/sec) 


2196633.860 

5643.485 

3361.685 

1410.426 


1000 meters 
1000 meters 
1000 m/sec 
1000 m/sec 
1000 m/sec 


At epoch, the satellite topographic coordinates are: 


Latitude = 

Longitude = 

Altitude = 

Velocity = 

Velocity azimuth = 
Flight path angle = 


16. 673 degrees North 
61. 501 degrees West 
1325537. 562 meters 
6718. 568 m/sec 
77. 739 degrees 
1. 291 degrees 


The trace of this trajectory over the Caribbean is shown in Figure 1 and the 
elevation angles vs time for each station are plotted in Figure 2. Note that the 
pass is almost overhead for Antigua. For this reason, Antigua is always chosen 
for the one station cases. 


3. MEASUREMENTS 


Data are simulated for each laser at 40 second intervals and for the altime- 
ter at 50 second intervals for Cases 1, 2, and 4a; at 10 second intervals for the 
lasers and the altimeter for Case 3; and at 1 second intervals for the lasers for 
Case 4b. These laser data rates fall between the one per second and one per 
minute data rates used by the Goddard and SAO lasers respectively. The cam- 
era data, right ascension and declination, are generated at 4 second intervals 
with a total time span of 24 seconds, thus simulating a GEOS 7-flash sequence. 
The camera data in all cases are equivalent to 1 plate per camera observed near 
the middle of the laser data s^an. The altimeter data rates are kept low to min- 
imize the influence of these data on the solution. 
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Simulated data are generated for the lasers, the cameras, and the altime- 
ter with the following RMS noise values: 


Laser range 
Laser angle 
Camera R.A. and Dec. 
Altimeter 


2 meters 
100 arc seconds 
1 arc second 
10 meters 


4. ERROR MODEL 


For each station that tracks , a range bias and the station survey are re- 
covered. An altitude bias is recovered for the altimeter. The camera and laser 
angles are assumed to be free of bias errors, although, for small quantities of 
data, the biases in these measurements may be as large as the RMS noise values 
assumed above. The a priori bias uncertainties for the error models are: 


Laser range bias = 2 

Altimeter range bias = 100 

Survey: Longitude, E =30 

Latitude, N =30 

Height (MSL) , V = 1 


meters 

meters 

meters 

meters 

meter 


IV. RESULTS 

All of the computer runs, except for Cases 3c and 4b, were made on the IBM 
System 360/91 computer using the NAP -II computer program operating in the 
simulation mode. This mode of operation utilizes the given orbital elements and 
measurement RMS noise values to generate simulated observations, forms the 
observation residuals , and outputs the error model recovery statistics for a 
given tracking situation. 

In ^il of the simulations , recovery is made for the six orbital elements , the 
survey for each station, a range b^as for each laser, and the altimeter bias. 

Initially , it was thought it might he necessary to tightly constrain either the 
orbital elements or the survey for at least one station in order to obtain a con- 
vergent solution. This question was investigated for the 4-laser, 0-camera 
(4LOC) case by obtaining solutions with various orbit and survey constraints (see 


7 



Table 1). The first soluaon used both a more constrained orbit (10m, 0.5m/s) 
and a fixed survey for one of the laser stations (GTK). The second solution held 
the same orbit constraints (10 m, 0. 5 m/s), but allowed for recovery of survey 
for all 4 stations. This increased the uncertainty in the altimeter bias recovery 
from 3.13 meters to 3. 28 meters. The third solution held the survey for one 
station fixed and relaxed the orbit constraints (1000m, 1000 m/s). This further 
increased the altimeter bias recovery uncertainty to 3.33 meters. Finally, the 
fourth solution used the relaxed orbit constraints (1000 m, 1000 m/s) and also 
recovered the survey for all 4 stations. This solution resulted in an altimeter 
bias uncertainty of 3. 42 meters. For these runs, the difference between the 
most constrained and the least constrained cases was only 0.29 meters in alti- 
meter bias recovery uncertainty. On this basis, it was decided that the a priori 
orbit and survey constraints could be relaxed with little penalty. 

Some additional constrained orbit solutions were run to evaluate the effect of 
adding cameras and the effect of laser timing uncertainties, (see Table 2). All 
of the solutions in Table 2 utilized the same orbit constraints (10m, 0. 5 m/s). 

The network configuration and constraints for the first two solutions in Table 
2 are identical to those for the first two solutions in Table 1, except that now a 
camera at GTK has been added to the 4 laser network. The first solution used 
a fixed survey for GTK. Adding a camera at GTK improved the altimeter bias 
recovery from 3. 13m (Table 1) to 2. 68 m (Table 2). For the second solution, 
allowing for recovery of survey at all 4 sites , adding a camera at GTK improved 
the altimeter bias recovery from 3.28m (Table 1) to 3.19m (Table 2), showing 
that adding a single camera to the 4 lasers is not as much help when the camera 
station survey is allowed to adjust. 

In solutions 3, 4, and 5 the number of cameras was increased to 2, 3 and 

4 respectively, with a corresponding reduction in altimeter bias recovery to 
3. 04m, and 2. 97 m, and 1. 54m. The addition of CUR camera data in solution 

5 was quite helpful, reducing the altimeter bias error to 1. 54m. This was prob- 
ably because of the earlier acquisition of data by CUR for the chosen trajectory. 

In solution 6, the configuration was the same as in solution 5, but now the 
lasers were allowed to have timing biases with a 1 ms uncertainty. This in- 
creased the altimeter bias recovery uncertainty to 2. 84m. For the remaining 
solutions, timing biases are not allowed, because station clocks can easily be 
synchronized in the Caribbean to within a few microseconds by a variety of tech- 
niques (e.g. , Loran-C). 

In summary, based on the results with these initial runs, the following as- 
sumptions are made in the solutions for the four cases outlined. The orbital 
elements are essentially unconstrained with an a priori uncertainty of 1000 m in 
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each position component and 1000 m/s in each velocity component. The station 
survey is recovered for all stations with a horizontal uncertainty in latitude and 
longitude of ±30 m for each, and a vertical uncertainty (relative to MSL) of ±1 m. 

This assumption for the vertical uncertainty is realistic only near the island 
tracking site, unless additional data are available. Away from the island, one 
must know the spatial mean sea level undulations relative to the station. These 
can be estimated from detailed ocean surface gravity measurements. Also, to 
some extent, one must know the sea state, the surface barometric pressure, and 
the surface winds. Many of these data can be made available for the Caribbean. 
However, the assumption of only a 1 -meter error in station height relative to 
the height of the spatial undulations in mean sea level may be too optimistic. 
More recent estimates, based on available ocean surface gravity data, quote 
this error as 2 or 3 meters. On the other hand, the laser range biases are as- 
signed an uncertainty of ±2m, which is probably too pessimistic, since more re- 
cent estimates, based on hardware improvements, quote a value of 0. 5m. Since 
these two parameters are highly correlated, the a priori estimate of one can be 
increased by 1.5 m and the other decreased by 1.5 m with probably little effect 
on the results for the altimeter bias. The altimeter bias is loosely constrained 
with an uncertainty of ±100 m. 

One might anticipate that the altimeter bias recovery from these solutions 
will be best for the high elevation angle data and will amount to approximately 
2.2 m, the root sum square of the 2 m laser range bias uncertainty and the 1 m 
height uncertainty. 


CASE RESULTS 
Case 1 : 

The purpose of this case is to determine how the number of lasers and cam- 
eras affects the altimeter bias recovery. Table 3 is a tabulation of the re- 
sults obtained by varying the number of lasers from a minimum of one to a 
maximum of four while varying the number of cameras from zero to the 
number of lasers, as described earlier. From Table 3, the estimated er- 
rors in recovery of altimeter bias , orbital elements , and station latitudes 
and longitudes are plotted in Figures 3, 4, 5, 6, and 7. 

From Figure 3 it appears that very little accuracy in altimeter bias recov- 
ery is gained beyond the 2 lasers, 2 cameras (2L2C) configuration. This is 
also evident in the recovery of the orbital elements , shown in Figures 4 and 
5 and in the horizontal survey recoveries shown in Figures 6 and 7. Dou- 
bling o f the tracking equipment from the 2L2C to the 4L4C network gains 
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only about 0. 3 meter accuracy in determination of altimeter bias. In addi- 
tion, it is highly unlikely that all four stations in the larger network will 
obtain data on the .same 1 pass, due to the weather dependence. 

The simplest and cheapest tracking station configuration which meets the 
4.1m requirement, is the 1L1C combination. Next simplest are the 2L1C, 
2L2C, and 3LOC configurations. 

However, as indicated in Reference 11, the most effective configuration is 
not necessarily the simplest. Effectiveness involves additional factors such 
as the probability of obtaining data and the duration of the requisite accu- 
racy in altitude determination. 

Estimates of the effect of weather and daylight on the relative probabilities 
of collecting acceptable altimeter calibration data with the various laser and 
camera configurations are developed in the Appendix. The results indicate 
that the 2L2C configuration is the most cost effective configuration having 
a reasonable probability of obtaining data. 

Case 2 : 

In case 2, it is assumed, for most of the runs, that lasers are available 
without cameras , and that laser angles provide the only directional infor- 
mation. The laser angle data is assumed to have an uncertainty of 100 arc 
seconds in azimuth and elevation. Simulations are made for 1, 2, 3, and 4 
lasers tracking without cameras, and also for the 4L1C configuration. The 
results are tabulated in Table 4 and the altimeter bias errors are plotted 
along with the Case 1 results in Figure 3. 

The one station 1LOC case (Antigua) is greatly improved by using the laser 
angles rather than the laser ranges alone. The altimeter bias uncertainty 
improves from ±38. 7 meters to ±4.7 meters with the addition of the laser 
angles. The importance of the laser angle data in this case due primar- 
ily to the high 87° maximum elevation angle of the pass over Antigua. When 
the one station situation is repeated using the Grand Turk station with laser . 
range and angles , and a maximum elevation angle for the pass of only 54° , 
the altimeter bias recovery ' ±25. 5m. The 2LOC case is improved from 
5. 1 m to 4. 4 m by adding laser angle data. 

As can be seen from Figure 3, when more than 2 lasers are tracking, the 
results are not changed by the addition of laser angle data. For only one or 
two lasers and no cameras, the laser angles are some help, but with cam- 
era angles the improvement in altimeter bias recovery is more pronounced. 
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Case 3: 


a In Case 3a, recovery of the altimeter bias is simulated using laser and 
camera data from Antigua at various elevation angles. For each of the 14 
one -minute mini -arcs between the minimum elevation of 16° and the maxi- 
mum elevation of 87° over the Antigua station, the state vector, altimeter 
bias, laser range bias, and station survey are recovered. The recovery 
uncertainties are given in Table 5. As anticipated, the uncertainty in the 
altimeter bias decreases as the elevation angle at the station increases. 

This is illustrated in Figure 8. The asymmetry between the ascending and 
descending mini -arc altimeter bias recoveries in Figure 8 is due to plotting 
the bias recoveries at the elevation angle of the beginning of each mini -arc, 
rather than at a more representative elevation angle near the middle of each 
mini -arc. The best mini -arc altimeter bias recovery for the 1L1C case is 
worse than the 1L1C altimeter bias recovery over the whole arc in Case 1 
because there are fewer observations in the mini -arc with which to deter- 
mine the 11 unknown parameters. 

b. In Case 3b, the objective is to determine the dependence of the altimeter 
bias recovery on elevation for combinations of stations. The computer sim- 
ulations vary the equipment from the maximum 4L4C case to the minimum 
1LOC case with several intermediate cases. An independent altimeter bias 
is recovered for most of the 14 mini-arcs. The orbit, survey, and bias 
parameters are recovered over the whole pass. This, in effect, gives the 
uncertainty of the altimeter bias recovery as a function of the network ele- 
vation angle and represents the situation when only short spans of altimeter 
data are obtained. 

The results are given in Table 6 and Figure 9. The five curves in Figure 9 
represent the altimeter bias uncertainty as a function of time. The elevation 
angles at each station as a function of time are given in Figure 2. Indicated 
on Figure 9 are the times of maximum elevation for the stations. These re- 
sults show the slight advantage for the 1L1C case of recovering most of the 
parameters over the whole pass rather than over each mini -arc as in Case 
3a. The plots also indicate that the larger the network and the better the 
tracking geometry the better the altimeter bias recovery at lower elevation 
angles. In the area of maximum elevation, the 4-station network provides 
only 0. 5 meter better accuracy than the 1 -station alone, provided the laser 
angles at the single station are used. Without the angle data the single sta- 
tion results are much worse, as can be seen in Table 6. 

c. In Case 3c, a simpler approach is employed. For the same trajectory 
points used in Case 3a, the altitude error is determined through a geomet- 
ric error propagation, assuming a survey uncertainty at ANT of (E, N, V) = 
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(30, 30, l)m, a range measurement uncertainty of 2 m, and an angle meas- 
urement uncertainty of either 1, 2, 4, 10, or 100 arc -sec. The results are 
shown in Figure 10. Again, it is evident that the altitude is determined best 
at the highest elevation angles and at these angles is not very dependent on 
the angular accuracy. For the lower elevation angles the altitude is better 
determined with the better angles. 

Case 4 : 

a. In Case 4a, the trajectory is tracked by 3 lasers at Antigua, Key West, 
and Panama. All recovered parameters and constraints are the same as 
in Case 1 (Table 3), where the same trajectory is tracked by the 3 more 
closely spaced lasers at Antigua, Grand Turk, and Curacao. 

The results are given in Table 7 , along with the earlier results from Table 
3 for comparison. The ANT, KWT, PAN network achieves an altimeter 
bias recovery of 4. 13m compared to 3. 80m with the ANT, GTK, CUR net- 
work. Thus, for this set of inputs, the simulations achieve a 0. 33m better 
altimeter bias recovery with the more closely spaced network. 

Both network results are dependent on the inputs and constraints. With a 
more constrained orbit (10 m, 0.5m/s), the ANT, KWT, PAN network achieves 
a 3. 33 m altimeter bias recovery. With a more constrained horizontal sur- 
vey of, say, 5m rather than 30m, this network would undoubtedly improve 
even further, but this run was not attempted. 

Both network results are also dependent on the high elevation angle data 
from Antigua for the chosen trajectory. For example, if TRD replaces ANT 
in the closely spaced network, the highest elevation angle data now comes 
from CUR at 74° , and the altimeter bias recovery , even for the constrained 
orbit case (10m, 0.5 m/s), deteriorates to 5.6m. 

b. In Case 4b, the ORAN orbital error propagation program is used to iden- 
tify the critical error sources in the determination of satellite height with 
simultaneous data from 2 lasers and no cameras (2LOC), and with more op- 
timistic constraints. 

The same orbital elements are used, but the epoch is changed to bring the 
trajectory approximately midway between the 2 lasers at ANT and TRD, 
providing 18 minutes of simultaneous tracking above 5° . The laser data 
rate is increased to one per second and the rms noise and bias a priori 
sigmas are both decreased to 1 m. The TRD survey is held fixed, and the 
ANT a priori survey error is set at (E, N, V) = (6, 6, 2)m, rather than 
the previous (30, 30, l)m. In this case the gravity field errors are also 
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modelled, using an a priori error estimate of 1 part in 10 3 * * 6 for GM, and the 
full difference (except for GM) between the SAO-M1 and the SAO COSPAR 
coefficients in the spherical harmonics expansion for the gravity potential. 

The results near the beginning, middle, and end of the pass are given in 
Table 8. Again, the satellite altitude determination is best near the middle 
of the pass, where elevation angles are highest. At this point the error in 
height is 1. 64 m, for the more optimistic constraints employed in this error 
propagation. With the less optimistic constraints used earlier, this com- 
pares with an altitude error of 2.2 m for the 1L1C error propagation in Case 
3c. No error propagations were done for the 2LOC configuration using the 
earlier constraints. 

When altimeter bias recovery, rather than error propagation, is attempted 
with the less optimistic constraints as in Cases 1,2, and 3a, the 2LOC con- 
figuration yields a bias recovery of 5. lm, 4.4m, and 6.2m respectively. 

To better evaluate the significance of the improvement in altitude determina- 
tion indicated in the 2LOC ORAN error propagation, the same more optimis- 
tic constraints should be applied consistently to error propagations and bias 
recoveries for the other laser and camera network configurations. 


V. CONCLUSIONS 

1. Considering the results of all the cases described, the simplest config- 
uration adequate for the altimeter bias determination appears to be the 1L1C 
combination. Assuming a 50% chance of clear weather, as needed by both the 
laser and the camera, and a 50% chance of darkness, as needed by the camera, 
the chance of obtaining data with this combination on random passes is about 25%. 
A drawback with this configuration, however, is that it provides sufficiently ac- 
curate altimeter calibration data only for high elevation passes above about 60° , 
and only for short durations, up to about 4 minutes. 

2. Adding more lasers and cameras improves the accuracy of the bias de- 
termination, extends the accuracy over a longer part of the pass and to lower 
elevation angles , and increases the chances of obtaining data. 

3. The 2LOC configuration, with the superior lasers and accurate relative 

station surveys indicated in Case 4b, is the next simplest configuration which is 

adequate. But, for realistic horizontal survey uncertainties, this configuration 

is probably less accurate than the 1L1C configuration for a high elevation pass. 
Assuming as before, a 50% chance of clear weather at either station, with no 
intersite weather correlation, and with no dependence on darkness, the chance 
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of obtaining data simultaneously from these stations is about 25%, the same as 
for the 1L1C case. 

4. By adding a camera to each of the 2 laser stations, we achieve a 2L2C 
configuration, which includes the acceptably accurate 2LOC and two 1L1C con- 
figurations. With the same 50% assumptions regarding the weather and dark- 
ness, the chance of obtaining data from the 2LOC subset is 25% as before, with 
half of these passes in darkness and suitable for data collection by the more 
accurate 2L2C network. In addition, there are cases when either one or the 
other 1L1C station can observe, increasing the overall probability of useful data 
collection to 50%, as shown in the Appendix. On the basis of the results of this 
simulation, the 2L2C configuration is probably the most effective combination 
from the standpoint of cost, accuracy, and probability of useful data acquisition. 


VI. RECOMMENDATIONS 

1. A network of two stations, each having a laser and a collocated camera 
(2L2C configuration), and situated for simultaneous observations (about 500 to 
1000 kilometers apart) is recommended as a basis for a satellite altimeter cali- 
bration such as that anticipated on GEOS-C. 

2. Some of these cases should be repeated using more optimistic horizontal 
survey and other constraints such as were used in Case 4b. This would benefit 
the multiple station solutions somewhat. 

3. The effect of using Doppler data in the short arc solutions should be in- 
vestigated. Such data will be available on GEOS-C in the Caribbean Sea area. 

4. The effect of using data from other ground stations on longer arcs, as 
described in reference 12, should be further investigated using a more optimistic 
assumption for the gravity field errors. Also the separately examined range, 
range rate, and angle data should be tried in various combinations. 

5. The effect of using satellite -to-satellite tracking, as described in refer- 
ences 13 and 14, should be further investigated using a realistic assumption for 
the gravity field errors. This technique seems to be potentially very accurate, 
provided the orbital arc extends over at least one revolution. This technique, 
because of its coverage where there are no ground stations, is valuable in con- 
junction with the altimeter for monitoring the altimeter spacecraft motions and 
thereby helping to separate variations in the altimeter data due to spacecraft 
motion from those due to ocean topography. However, since the satellite-to- 
satellite technique has not yet been evaluated using real data, it would be risky 
to rely on this technique alone for the GEOS-C satellite altimeter calibration. 
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APPENDIX 


PURPOSE: 


To estimate the effect of weathe r and daylight on the probability of useful 
data collection by various laser and camera configurations for altimeter cali- 
bration purposes. 

Assumptions : 

1. All configurations more complex than 1LOC can obtain useful data for alti- 
meter calibration purposes (based on Case 1 and 4b results). 

2. Satellite passes occur with equal probability over each station. If a pass is 
above the elevation horizon for one station, it is above this horizon for the others. 

3. Daylight (L) or night (N) passes occur with equal probability P(L) = : P(N) = 
0.5. If a pass is light or night for one station, it is the same for the others. 

4. The weather correlation between stations is assumed to be zero. Transpar- 
ent (T) or Opaque (O) passes occur with equal probability, P(T) = P(O) = 0. 5. 

A second case with probability P(T) = 0.4, P(O) = 0. 6, which is more repre- 
sentative of the GEOS-1, 2 results with MOTS cameras in the Caribbean (ref- 
erence 15) is also examined. 


DATA COLLECTION PROBABILITIES 


For the equally likely probabilities chosen for transparent or opaque skies, 
the desired data collection probabilities for the various station configurations 
can be derived as the outcome of coin tossing experiments , where the coins are 
labelled transparent (T) or opaque (O), rather than heads or tails. Furthermore, 
when cameras are involved and darkness is a factor, the equally likely probabil- 
ities chosen for darkness or daylight can be represented as another coin in the 
experiment labelled night (N) or light (L). 

One Station 


Network A 0 = 1LOC (Data not adequate) 

This network requires transparent sky passes to acquire data. The equally 
likely coin toss outcomes are T or O. Thus P(A 0 ) = P(T) = 1/2. 
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Network A! - 1L1C (Simplest case with adequate data) 

This network requires both night passes and transparent sky passes in or- 
der to acquire data. The equally likely outcomes, when tossing two coins, one 
labelled N or L, the other T or O are,' 

N T 
N O 
L T 
L O 

The desired result (N, T) has a probability P(A t ) = P(N, T) = 1/4 - 
Two Stations 

Network B 0 ~ 2LOC 

This network requires transparent sky passes at both laser sites simulta- 
neously, but is independent of whether the passes are day or night. The equally 
likely outcomes when tossing two coins, both labelled T or O are: 

T T 
T O 
O T 
O O 

P(B 0 ) = P(T, T) = 1/4 
Network Bi = 2L1C 

This network includes both the 2LOC and the 1L1C cases considered above. 
Considering the night and light passes separately, the possible weather outcomes 
at the two sites are: 

Night (2 L1C or 1L1C) Light (2LOC) 

T T T T 

TO TO 

O T O T 

O O O O 
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At night, acceptable data can be obtained with the 2L1C systems if it is clear at 
both sites, or with the 1L1C systems if it is clear at the site with the camera. 

It is implicitly assumed that at night a possible 2LOC pass is equivalent to a 
2L1C pass, since both cases require the same clear skies. P(B, /N) = PfT, T) 
+ P(T, O) - 1/4 + 1/4 = 1/2 

In the day passes, acceptable data can be obtained with the 2LOC systems 
if it is clear at both sites. P(Bi /L) = P(T, T) = 1/4 

The total probability of obtaining data with this network is: 

P(B i ) = P(B i/N)P(N) + P(Bi/L)P(L) = (1/2) (1/2) + (1/4) (1/2) = 3/8 

Network B 2 = 2L2C 

Since this is still a two station case, the possible weather outcomes given 
above apply. On the night passes, acceptable data can be obtained with the 2L2C 
systems if it is clear at both sites, or with either of the 1L1C systems, if it is 
clear at either site. P(B 2 /N) = P(T, T) + P(T, O) + P(0, T) = 3/4 

P(B 2 /L) = 1/4 as before 

P(B 2 ) = P(B 2 /N)P(N) + P(B 2 /L)P(L) = 3/4 • 1/2 + 1/4 • 1/2 = 1/2 
Three Stations 


Proceeding as above the day or night possible weather outcomes are: 
T TT 
T TO 
T OT 
T OO 
O TT 
O TO 
O OT 
O OO 
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Network C ( , - 3LOC 


This requires clear weather at all three sites or at any two (2LOC) sites for 
adequate data, day or night. 

P(Co) - P(T,T,T) + P(T,T,()) + P(T,O.T) + P(0,T,T) = 1/8 + 1/8 + 1/8 + 1/8 = 1/2 
Network C t - 3U C 

The night passes can provide data for the 3L1C, 2L1C, 2LOC, or 1L1C sys- 
tems and therefore require clear weather at all three sites , or at any pair of 
sites, or at the camera site alone. 

P(Ci/N) - P(C 0 /N) + P(T ,0,0) = 1/2 + 1/8 = 5/8 

The day passes provide data for the 3LOC or 2LOC systems as above. 

P(Ci/L) = P(C 0 ) - 1/2 

P(Ci) - P(C!/N)P(N) + P(C | /L)P(L) = 5/8 • 1/2 + 1/2 • 1/2 = 9/16 
Network C 2 - 3L2C 
Proceeding as above 

P(C 2 /N) = P(Ci/N) + P(0,T,0) = 5/8 + 1/8 = 3/4 
P(C 2 /L) = 1/2 again 

P(C 2 ) = P(C 2 /N)P(N) + P(C 2 /L)P(L) = 3/4 • 1/2 + 1/2 • 1/2 = 5/8 
Network C 3 = 3L3C 

P(C 3 /N) = P(C 2 /N) + P(0,0,T) = 3/4 + 1/8 = 7/8 
P(C 3 /L) = 1/2 

P(C 3 ) = P(C 3 /N)P(N) + P(C 3 /L)P(L) = 7/8 • 1/2 + 1/2 • 1/2 = 11/16 

1 

Four Stations 

Network D 0 = 4LOC 
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This requires clear weather for any two of the four sites. There are 16 
possible weather events, the 8 previously given for 3 stations when the fourth 
station is clear, and these 8 again when the fourth station is cloudy. 

P(D n ) = 7/16 + 4/16 = 11/16 

Network Di = 4L1C 

P(D,/N) = P(D 0 /N) + P(T, 0,0,0) 

= 11/16 + 1/16 = 3/4 

P(Dj/L) = P(D 0 /L) = 11/16 

P(Di) = P(Di/N) P(N) + P(D!/L)P(L) = 3/4 • 1/2 + 11/16 . 1/2 = 

23/32 = P(Do) + 1/32 

Network D 2 = 4L2C 

P(D 2 /N) = P^/N) + P(0,T,0,0) = 3/4 + 1/16 = 13/16 
P(D 2 L) = P(D 0 /L) = 11/16 

P(D 2 ) = 13/16 • 1/2 + 11/16 • 1/2 = 3/4 = P(D t ) + 1/32 

Network D3 = 4L1C 

P(D 3 ) = P(D 2 ) + 1/32 = 25/32 

Network D 4 = 4L4C 

P(D 4 ) = P(D 3 ) + 1/32 = 13/16 

The above data probabilities were derived assuming a probability of clear 
weather P(T) = 50%. By slightly more involved techniques the data probabili- 
ties can be derived assuming P(T) = 40%, which is more representative of the 
probabilities actually experienced with MOTS cameras tracking GEOS-1, 2 from 
Jamaica and Puerto Rico. The results for both weather probabilities are sum- 
marized on the next page. 

Note that for P(T) = 50%, the probabilities of obtaining useful data P(Ai) = 

P(B 0 ); P(B 2 ) = P(C 0 ); P(C 3 ) = P(D 0 ). Each time another laser is added to the 

network, it serves to replace all the cameras in the previous network with re- 1 

spect to data collection capability as defined in this report. Thus, the 25% data 
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Network 

Configuration 

Useful Data Probabilities 

P(T) = 50% 

P(T) = 40% 

Ai 

1 Station 

1L1C 

1/4 = .250 

.200 

B 0 

2 Stations 

2LOC 

2/8 = .250 


Bi 


2L1C 

3/8 = .375 

.280 

b 2 


2L2C 

4/8 = . 500 

.400 

Co 

3 Stations 

3LOC 

8/16 = .5000 

.352 

c, 


3L1C 

9/16 = .5625 

.424 

C 2 


3L2C 

10/16 = .6250 

.496 

c 3 


3L3C 

11/16 = .6875 

.568 

Do 

4 Stations 

4LOC 

22/32 = .68750 

. 5248 

D i 


4L1C 

23/32 = .71875 

.5780 

d 2 


4L2C 

24/32 = .75000 

.6312 

d 3 


4L3C 

25/32 = .78125 

.6844 

d 4 


4L4C 

26/32 = .81250 

.7376 


collection probability of a 2LOC configuration can be doubled by adding 2 cam- 
eras or 1 laser. Similarly, the 50% probability with a 3LOC configuration can 
be increased to 68. 75% by adding 3 cameras or 1 laser. 

For the probabilities with P(T) = 40%, the data collection probabilities are, 
of course, lower, but the cameras are relatively more effective in improving 
these probabilities. 

Assuming the altimeter satellite carries both laser corner reflectors and 
optical beacons , the 2L2C configuration is probably the most economical ground 
station configuration having a reasonably large probability for collection of suit- 
ably accurate altimeter calibration data and is therefore recommended as the 
most effective laser and camera network for this purpose. 
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Table 1 


Exploratory Results with Different Survey and Orbit Constraints 


Solution 

1 

2 

3 

4 

Configuration 

4L0C 

4L0C 

4L0C 

4L0C 

A priori orbit (m, m/s) 

m 

o 

o 

10, 0.5 

1000, 1000 

1000, 1000 

Survey (Fix, Find) 

1, 3 

0, 4 

1. 3 

0. A 

Alt. Bias H(m) 

3.13 

3.28 

3.33 

3.42 

Orbit (m, m/s) 





E, E 

1.7, .020 

4.3, .04 

2.1, .04 

9.9, .05 

F , t 

5.7, .003 

8.4, .01 

10.7, .006 

17.5, .02 

G, G 

8.1. .041 

9.0, .01 

18.4, .08 

25.1. .15 


Range Bias. B(m) & Vert. Survey V (m) 


ANT 

(B, 

V) 

1. 

5, 

.97 

1.6, 

•9* 

1.5, 

.98 

1. 

6, 

.99 

GTK 

(B, 

V) 

1. 

■ 8, 

- 

1.9, 

.99 

1.9, 

- 

1. 

9, 

.99 

CUR 

(B, 

V) 

1. 

7, 

.99 

1.7, 

.99 

1.7, 

.99 

1. 

8, 

.99 

TRD 

(B, 

V) 

1. 

,8, 

.99 

1.8, 

.99 

1.8. 

.99 

1. 

8*. 

.99 

Horiz. 

Survey (m) 











ANT 

(E, 

N) 

7, 

.2, 

27.3 

9.4, 

27.9 

11.3, 

27.6 

16. 

o, 

28.4 

GTK 

(E, 

N) 

- 


- 

12.3, 

19.5 

- 

- 

22. 

o, 

20.9 

CUR 

(E, 

N) 

8, 

.3, 

13.1 

11.2, 

16.6 

14.7, 

15.3 

16. 

9, 

20.2 

TRD 


JD 

8, 

‘Az, 

9.7 

13.5, 

10.0 

17.1, 

17.9 

19. 

o*_ 

23.1 


A priori parameter uncertainties for all solutions : 

Survey: (E, N, V) = +(30, 30, l)m 

Laser Range Bias : B = +2m 

Altimeter Bias : H = +100m 
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Table 2 


Exploratory Results with Different Camera Configurations 


Solution 

1 

2 

3 

4 5 

6 

Configuration 

4L1C 

4L1C 

4L2C 

4L3C 4L4C 

4L4C 

Camera sites 

GTK 

GTK 

GTK, ANT or TRD 

GTK, ANT, TRD 


A priori Orbit 

10, 0.5 

10, 0.5 

10, 0.5 

10, 0.5 10, ■ ■ . 5 

10, 0.5 

(m, m/s) 






Survey (Fix, Find) 1,3 

0,4 

0,4 

0,4 0,4 

0,4 

A priori Laser 

Time Bias 

(ms 1 , ... 



1.0 

Alt. Bias, H(m) 

2.68 

3.19 

3.04 

2.97 1.54 

2.84 

Orbit (m, m/s) 




■ 


E, E 

1.7, .008 

3.5, .080 

1.8, .006 

1.7, .005 0.6, .002 

1.9, .004 

F, F 

2.2, .003 

4.8, .004 

2.4, .003 

2.0, .003 0.9, .001 

2.0, .002 

G, G 

5.8, .016 

8.8, .027 

3.5, .016 

2.3, .014 1.1. .003 

2.3. .010 


Range Bias, B(m), and Vert. Survey, V (m) , Time Bias (run 6 only), T(ms) 


ANT 

GTK 

(B ,V) 
(B,V) 

1.4, 

1.6 

.96 

1.6, 

1.9, 

.98 

.99 

1.5, 

1.9, 

.97 

.99 

1.5, 

1.9, 

.97 

.99 

1.3, 

1.9, 

.94 

.99 

T(ms) 
1.5, .96, .98 
1.9, .99, .98 

CUR 

(B,V) 

1.7, 

.99 

1.7, 

.99 

1.7, 

.99 

1.7, 

.99 

1,7, 

.99 

1.7, .99, .98 

TRD 

(B ,V) 

1.8, 

.99 

1.8, 

.99 

1.8, 

.99 

1.8, 

.99 

1.8. 

.99 

1.8. .99. .98 

Horiz. 

Survey 

(m) 











ANT 

(E,N) 

6.1, 

26.2 

6.4, 

27.7 

6.0, 

27.0 

5.9, 

26.9 

5.5, 

24.8 

8.0, 26.5 

GTK 

(E,N) 

- 

- 

5.7, 

6.8 

2.9, 

5.3 

2.6, 

5.3 

1.7, 

4.0 

5.8, 5.0 

CUR 

(E,N) 

3.6, 

17.0 

4.9, 

11.8 

3.5, 

11,0 

3.3, 

10.6 

2.5, 

8.5 

6.2, 10.4 

TRD 

Mi N). - 

3.2, 

7.8 

4.2, 

9.6 

_2jJLt_ 

6.8 

2.4> 

6.2 

1.6, 

4.8 

5.7. 6.1 


A priori parameter uncertainties for all solutions: 

Survey: (E, N, V) = +(30, 30, l)m 

Laser Range Bias : B = +2m 

Altimeter Bias : H = +100m 

Time Bias (run 6 only) : T = +lms 
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Table 3 


Case 1 


Lasers 

4 

4 

4 

4 

4 

3 

3 

3 

3 

2 

2 

2 

1 

1 

Cameras 

4 

3 

2 

1 

0 

3 

2 

1 

0 

2 

1 

0 

1 

0 

Alt. Bias 














1m) 

2.87 

2.99 

3.14 

3.36 

3.42 

3.00 

3.16 

3.40 

3.80 

3.20 

3.50 

5.10 

3.60 

38.70 

Orbit (: 

m, m/s) 














E 

1.6 

1.9 

2.1 

2.8 

9.9 

1.9 

2.1 

2.9 

12.9 

2.1 

3.4 

17.6 

4.6 

250. 

F 

2.0 

2.6 

2.7 

2.9 

17.5 

2.6 

2.8 

3.0 

20.1 

2.8 

3.0 

30.1 

3.4 

90.4 

G 

2.3 

2.8 

2.9 

6.7 

25.1 

2.8 

2.9 

7.4 

42.5 

2.9 

8.1 

76.2 

13.7 

827. 

E 

.004 

.005 

.006 

.044 

.054 

.005 

.006 

.049 

.085 

.006 

.055 

.144 

.097 

.752 

F 

.003 

.003 

.003 

.013 

.017 

.003 

.003 

.014 

.022 

.003 

.016 

.028 

.028 

.216 

6 

.010 

.011 

.015 

.124 

.149 

.011 

.015 

.138 

.220 

.015 

.157 

.334 

.276 

2.09 

Range Bias (m) 



* 











ANT 

1.5 

1.5 

1.5 

1.6 

1.6 

1.5 

1.6 

1.6 

1.6 

1.6 

1.6 

1.7 

1.7 

1.9 

GTK 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 


- 

CUR 

1.7 

1.7 

1.7 

1.7 

1.8 

1.7 

1.7 

1.7 

1.8 

- 

- 

- 

- 

- 

TRD 

1.8 

1.8 

1.8 

1.8 

1.8 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Survey 

(m) 














ANT E 

5.9 

5.9 

6.1 

6.4 

16.0 

5.9 

6.1 

6.4 

18.5 

6.1 

6.6 

28.9 

6.8 

29.9 

N 

26.6 

27.0 

27.4 

28.1 

28.4 

27.1 

27.5 

28.3 

28.5 

27.8 

28.6 

28.9 

28.9 

29.9 

V 

.97 

.97 

.98 

.99 

.99 

.97 

.98 

.99 

.99 

.98 

.99 

.99 

1.0 

1.0 

GTK E 

2.4 

2.7 

3.1 

11.5 

22.0 

2.7 

3.1 

12.7 

25.9 

3.1 

14.4 

28.8 

- 

_ 

N 

4.9 

5.2 

5.7 

16.8 

20.9 

5.3 

5.7 

18.5 

21.9 

5.8 

20,2 

29.4 

- 

- 

V 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

- 

- 

CUR E 

3.3 

3.7 

3.7 

8.8 

16.9 

3.8 

3.8 

9.7 

23.4 

• 


• 

• 

_ 

N 

10.4 

10.5 

10.9 

13.7 

20.2 

10.6 

10.9 

14.3 

25.2 

- 

- 

- 

- 

- 

V 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

.99 

- 

- 

- 

- 

- 

TRD E 

2.4 

2.8 

2.8 

11.7 

19.0 


_ 

. 

_ 

- 

- 

- 

_ 

- 

N 

6.0 

6.1 

6.3 

8.8 

23.1 

- 

- 

- 

- 

- 

- 

- 

- 

- 

V 

.99 

.99 

.99 

.99 

.99 

- 

- 

- 

- 

- 

- 

- 

- 

- 
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Table 4 


Case 2 


Lasers* 

Cameras 

4 

0 

3 

0 

2 

0 

I 1 

0 

i 2 

0 

4 

1 

Alt. Bias (m) 

3.40 

3.66 

4.37 

4.70 

25.5 

3.35 

Orbit (m, m/s) 

E 

8.7 

10.6 

12.6 

18.7 

18.4 

2.8 

v 

15.9 

18.3 

25.8 

28.0 

99.7 

2.9 

G 

23.2 

36.2 

55.3 

70.1 

211.9 

6.5 

t 

.052 

.077 

.114 

.21 

.43 

.043 

i 

.015 

.019 

.022 

.05 

.06 

.013 

G 

.140 

.194 

.264 

.56 

.88 

.121 

Range Bias (m) 

ANT 

1.6 

1.6 

1.7 

1.89 

- 

1.6 

GTK 

1.9 

1.9 

1.9 

- 

1 r 99 

1.9 

CUR 

1.8 

1.8 

- 

- 

- 

1.7 

TRD 

1.8 

- 

- 

- 

- 

1.8 


ANT E 

14.9 

17.2 

24.7 

26.7 

- 

6.4 

N 

28.4 

28.4 

28.8 

29.7 

- 

28.1 

V 

.99 

.99 

.99 

.99 

- 

.99 

GTK E 

19.6 

21.8 

22.7 

- 

29.8 

11.2 

•. N 

20.4 

21.4 

27.4 

- 

29.7 

16.4 

V 

.99 

.99 

.99 

- 

.99 

.99 

CUR E 

16.0 

21.8 

— 

- 

- 

8.6 

N 

19.6 

23.3 

- 

- 

- 

13.6 

V 

.99 

.99 

- 

- 

- 

.99 

TRD E 

18.2 

- 

_ 

- 

- 

11.4 

N 

21.4 

- 

- 

- 

- 

8.7 

V 

.99 

- 

- 

- 

- 

.99 


*Laser Range, Azimuth, and Elevation 
^■Antigua was the single station (best case) 
^Grand Turk was the single station (worst case) 
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Table 6 


Case 3b 


Lasers 

4 

3 

2 

* \ 

1* 

1 

Cameras 

4 

0 

L 

l 

0 

0 


Alt. Bias (m) 


Arc 1 

5.06 

6.83 

7.68 

8.92 

- 

- 

Arc 2 

- 

- 

- 

- 

8.18 

33.62 

Arc 3 

4.72 

5.97 

6.42 

7.16 

- 

- 

Arc 4 

- 

- 

- 

- 

6.54 

30.44 

Arc 5 

4.52 

5.33 

5.43 

5.74 

- 

- 

Arc 6 

4.49 

5.12 

5.09 

5.26 

5.39 

29.67 

Arc 7 

4.50 

4.98 

4.87 

5.01 

- 

- 

Arc 8 

4.55 

4.93 

4.79 

5.03 

5.14 

29.20 

Arc 9 

4 . 65 

4.96 

4.86 

5.29 

- 

- 

Arc 10 

4.78 

5.07 

5.04 

5.78 

5.89 

28.55 

Arc 11 

- 

- 

- 

- 

- 

- 

Arc 12 

5.16 

5.49 

5.72 

m 

OJ 

r-- 

7.34 

29.37 

Arc 13 

- 

- 

- 

- 

- 

- 

Arc 14 

5.64 

6.10 

6.65 

9.00 _ 

9.11 

34.80 

Orbit (m, m/s) 







E 

1.7 

13.6 

3.6 

5.3 

10,6 

195.7 

F 

1.9 

20.3 

2.8 

3.3 

20.8 

96.3 

G 

2.3 

44.5 

8.5 

14.6 

27.0 

65.9 

t 

.007 

.09 

.06 

.11 

.13 

2.7 

t 

.003 

.03 

.02 

.03 

.04 

0.8 

G 

.010 

.24 

.17 

.29 

.37 

7.3 

Ranee Bias (m) 







ANT 

1.2 

1.5 

1.7 

1.9 

1.9 

1.9 

GTK 

1.4 

1.5 

1.7 

- 

- 

- 

CUR 

1.3 

1.4 

- 

- 

- 

- 

TRD 

1.4 

- 

- 

- 

- 

- 

Survey (m) 







ANT E 

6.4 

18.5 

6.9 

7.8 

19.2 

29.9 

N 

27.5 

29.0 

29.1 

29.3 

29.5 

29.9 

V 

.97 

.99 

.99 

.99 

.99 

.99 

GTK E 

2.2 

27.2 

14.9 

- 

- 

- 

N 

4.3 

22.2 

20.7 

- 

- 

- 

V 

.99 

.99 

.99 

- 

- 

- 

CUR E 

3.4 

24.8 

- 

• 

- 

- 

N 

8.5 

25.0 

- 

- 

- 

- 

V 

.99 

.99 

- 

- 

- 

- 

TRD E 

2.3 

- 

- 

- 

- 

- 

N 

4.9 

- 

- 

- 

- 

- 

V 

.99 

- 

- 

- 

- 

- 


*Laser range, azimuth, elevation 
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Table 7 


Case 4a 


Configuration 

3L0C 

3LOC 

Laser Sites 

ANT, KWT , PAN 

ANT, GTK, CUR 

A priori orbit (m. m/s) 

1000m, lOOOm/s 

1000m, lOOOm/s 

Alt. Bias (m) 

4.13 

3.80 


Orbit: (m, m/s) 


E, E 


16.2, .062 


12.9, 

.085 

F, f 


20.4, .014 


20.1, 

.022 

G , u 


57.4, .143 


42.5, 

.220 

Range Bias , 

B(m) & Vert. 

Survey V (m) 




(B, V) 

ANT 

1.6, 0.99 

ANT 

1.6, 

.99 


KWT 

2.0, 1.00 

GTK 

1.9, 

.99 


PAN 

1.8, 1.00 

CUR 

1.8, 

.99 

(E, N) 

ANT 

28.7, 19.4 

ANT 

28.5, 

18.5 


KWT 

22.1, 27.9 

GTK 

21.9, 

25.9 


PAN 

24.3, 23.3 

CUR 

25.2, 

23.4 
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Error Contributions Near Beginning, Middle, and End of Pass Case 4b 
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TIME (min) 

Figure 2. Elevation Angles for 4-Station Network 



NOTE: DATA SIMULATED FOR LASERS AND CAMERAS WAS ORDERED BY STATION 
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Figure 4. Case 1, Satellite Position Recovery vs Equipment Tracking 
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Figure 7. Case 1, Satellite Longitude Recovery vs Equipment Tracking 
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Figure 8. Case 3a, Altimeter Bias Uncertainty vs Elevation Angle 
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Figure 9. Case 3b, Altimeter Bias Recovery vs Time of Pass 
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Figure 10. Case 3c, Altitude Uncertainty vs Elevation Angle 



